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We consider quasiparticle charge and energy imbalances in a thin superconductor weakly 
coupled with two normal-metal electrodes via tunnel junctions at low temperatures. Charge 
and energy imbalances, which can be created by injecting quasiparticles at one junction, 
induce excess tunneling current I cx at the other junction. We numerically obtain 

Jex as a 

function of the bias voltage Vdet across the detection junction. We show that I ox at the zero 
bias voltage is purely determined by the charge imbalance, while the energy imbalance causes 
a nontrivial Vdct-dependence of I cx . The obtained voltage-current characteristics qualitatively 
agree with the experimental result by R. Yagi [Phys. Rev. B 73 (2006) 134507]. 
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Quasiparticle current injected into a superconductor creates a nonequilibrium distribution 
of quasiparticles. 1-8 It has been shown that such a nonequilibrium distribution is decomposed 
into two components. 6 Let f(r, e) be the local deviation of a quasiparticle distribution from the 
equilibrium one, where e denotes quasiparticle energy measured from the chemical potential. 
One component is the number difference f<Y{r,e) = (/(r,e) — /(r,— e))/2 between electrons 
and holes, and the other is the average e) = (f(r, e) + f(r, — e))/2 of electron and hole 
numbers. The former one describes charge imbalance of quasiparticles, so we call it particle 
mode, and the latter one describes energy imbalance, so we call it energy mode. Alternatively, 
/t (/l) is called transverse (longitudinal) mode since it is coupled with phase (amplitude) 
variations of the pair potential. Note that charge imbalance can arise in a superconductor 
without breaking the charge neutrality since the corresponding variation of Cooper-pair den- 
sity compensates it. Charge imbalance induces excess quasiparticle current, which results in 
a potential difference between pairs and quasiparticles. 1, 3,4 Past experiments have detected 
charge imbalance by measuring this potential difference. 2 ' 5 In contrast, a clear experimental 
detection of energy imbalance has not been reported so far. 

Since the past experiments adopted the current balance method to detect charge imbal- 
ance, we could approach only the case where the bias voltage Vd e t of the detection junction is 
equal to zero. Recently, Yagi succeeded to measure the Vd e t-dependence of the tunneling cur- 
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rent through a superconductor-normal metal tunnel junction under quasiparticle injection at 
low temperatures. 8 The device employed in the experiment consists of a thin superconductor 
(Al) to which several thin normal-metal electrodes (Au) are connected via tunnel junctions. 
One junction is used as a quasiparticle injector, while the tunneling current is detected at 
another junction. One expects that charge imbalance created at the injection junction induces 
the excess tunneling current at the detection junction if the separation between the two junc- 
tions is much shorter than the corresponding relaxation length. Note that energy imbalance 
has been believed to be irrelevant to the excess current. If charge imbalance is absent, the 
tunneling current is given by 



where i?dct, A and /fd are the tunnel resistance of the detection junction, the energy gap and 
the Fermi-Dirac distribution function, respectively, and N± is the normalized BCS density of 
states given by 



excess current I ex arising from the charge imbalance. Naively speaking, J ex does not depends 
on Vdet unless | Vdet I A/e because its magnitude is determined by the total amount of the 
charge imbalance accumulated in the superconductor side of the junction. However, a notable 
deviation from /t(Vdet) = ^q(^det)+^ex is detected in the experiment although the result in the 
absence of quasiparticle injection is well explained by eq. (1). This indicates that the excess 
current has a nontrivial V^et -dependence, which cannot be explained by existing theories. 

In this letter, we study the excess tunneling current in a superconductor-normal metal 
junction under quasiparticle injection at low temperatures, to explain the experimental result 
by Yagi. We propose a simple model for the device employed in ref. 8, and numerically obtain 
the excess current function of the bias voltage Vdet taking account of not only the 

charge and energy imbalances in the superconductor but also nonequilibrium quasiparticles 
in the detection normal-metal electrode. We show that a nonequilibrium quasiparticle distri- 
bution created in the normal-metal electrode by the energy imbalance can contribute to the 
excess current although the energy imbalance itself has no direct contribution. We also show 
that / ex at the zero bias voltage is purely determined by the charge imbalance, while the energy 
imbalance indirectly causes a nontrivial V^et-dependence of I ex . The obtained voltage-current 
characteristics qualitatively agree with the experiment. We set h = &b = 1 throughout this 
letter. 

We consider a thin superconducting wire which is coupled with two thin normal-metal 
electrodes via tunnel junctions (see Fig. 1). Let Lg and Ln be the lengths of the superconductor 
and the normal-metal electrodes, respectively. We focus on the excess current through the right 
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Fig. 1. Schematic picture of the model system which consists of a thin superconductor (S) coupled 
with two thin normal-metal electrodes (N) via tunnel junctions. The lengths of the superconductor 
and normal-metal electrodes are Lg and Ln, respectively. The left and right junctions serve as 
quasiparticle injector and detector, respectively. 

junction in the presence of charge and energy imbalances created by quasiparticle current 
injection at the left junction. Let Vj e t an d Mnj be the bias voltages in the right and left 
junctions, respectively. In order to calculate the excess current, we need to obtain energy- 
dependent quasiparticle distributions in both the superconductor and the right normal-metal 
electrode for detection. A nonequilibrium quasiparticle distribution in the left normal-metal 
electrode for injection is not important for our argument, so we neglect it in the following. 
Assuming that the superconductor and the normal electrodes are very thin, we adopt a simple 
one-dimensional model for the device. We introduce the x axis in the superconductor on which 
the left and right junctions are at x = Xi n j and x = x^ ct , respectively, and the y axis in the right 
normal-metal electrode on which the right junction is at y = 0. In terms of the particle and 
energy modes, the nonequilibrium quasiparticle distribution in the superconductor is given 



where /fd(c) is the Fermi-Dirac distribution function with the chemical potential of the 
superconductor. We hereafter assume that the magnitude of the energy gap A is unaffected by 
quasiparticle injection everywhere in the superconductor. This allows us to consider /t,l(£, e) 
only for |e| > A. We turn to the nonequilibrium quasiparticle distribution in the right normal- 
metal electrode, which is perturbed by an electric field E and quasiparticle tunneling from 
the superconductor. If these perturbations are neglected, the quasiparticle distribution is 
simply given by g^(y,e) = /fd(c + eVdet)- Here and hereafter, we measure quasiparticle 
energy from in both the superconductor and the normal-metal electrodes. We consider 
roles of the perturbations. The field E induces a polarization of the distribution function in 
the momentum space. This polarization is not important for considering the tunneling current 
although it sustains the electron current in the normal-metal electrode. Thus, we neglect the 



by 6 



9s(x, e) = /fd(c) + h{x, e) + f L (x, e) 



(3) 
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effect of E. Accordingly, we approximately obtain 

9n(v, e) = /fd(£ + eV^ct) + / N (y, e), (4) 

where the second term represents the contribution from the quasiparticle tunneling. Exactly 
speaking, the charge neutrality in the vicinity of the junction is slightly broken under this 
assumption. To exactly ensure it, we need to introduce the shift of the chemical potential, 

My) = - J de /N(y,e), (5) 
into the first term. However, in actual situations, the shift is very small as long as the excess 
current is not so large. Furthermore, it does not affect the tunneling current unless |Vd c t| 
A/e. We have thus neglected 5^{y). As noted above, the tunneling current at the right junction 
is given by I t = I q + I cx . In terms of the distribution functions, we can express the excess 
current as 

4x = deN 1 (e)(2h(x dct ,e) - / N (0,e) - / N (0,-e)Y (6) 

where R^ et is the tunnel resistance of the right junction. Due to the asymmetric relation of 
e) = — e), 6 the energy mode does not appear in eq. (6). Note that the contribution 

from the quasiparticles in the normal-metal electrode is governed by /n(0, e) + /n(0, — e) and 
the excess current shows no Vdet-dependence if we completely neglect it. We shall show that 
/n(0, e) + /n(0, — e) is mainly determined by the tunneling of quasiparticles in the energy 
mode in collaboration with phonon-mediated energy relaxation. This means that the coupling 
between /n and the energy mode /l is essential to describe the Vdet-dependence. 

We present Boltzmann equations for /t, /l and /n, based on which we obtain /t and /n 
to calculate the excess current. The particle and energy modes in the superconductor obey 6 ' 9 

N?(e)D s ^fr(x, e) - -^Lf T ( x , e) + I T {x, e, {/t}) + Pr(x, e) = 0, (7) 

OX T conv {€) 

d 2 

D s -^f L (x,e)+I L (x,€,{f L }) +P L (x,e) = 0, (8) 

where D$ and r CO nv(e) represent the diffusion constant and the conversion time for charge 
imbalance, respectively. If impurity scattering with the pairing anisotropy is the dominant 
source, the conversion time is given by 10 

where ((ap) 2 ) ~ 0(0.01) is a parameter that characterizes the pairing anisotropy 11 and r the 
elastic relaxation time. The terms It and 1^, called collision integrals, describe influences of 
inelastic phonon scattering. They are given by 6 ' 9 

J T , L (s, e,{/}) =-2 J deV s (e,e , )M TiL (e,e / ) 

J cosh 2 f(x, e) - cosh 2 (j^) f(x, e') 



sinh(^j^) cosh(^) cosh(|^) 



+ 2f(x,e)f(x,e')\ (10) 
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with Mr(e,e') = N^N^e'), M L (e,e') = N^N^e') - R 2 (e)R 2 (e') and 

o"s(e, e') = ^psign(e' - e) • (e' - e) 2 , (11) 

where i?2(e) = (A/e)7Vi(e) and as characterizes the strength of the electron-phonon coupling 
in the superconductor. In eq. (10), the term with cosh 2 (e/2T)/(a;, e) (cosh 2 (e'/2T)/(x, e')) 
represents the so-called scattering-out (scattering-in) term. The nonlinear term with 
f(x, e)f(x, e'), which has been neglected in refs. 6 and 9, is important in obtaining the energy- 
mode /l at low temperatures and low energies, where inelastic phonon scattering is very weak. 
The coupling terms Pr and Pl represent the tunneling of quasiparticles between the super- 
conductor and the normal electrodes. For our purpose, we must consider the nonequilibrium 
quasiparticles induced in the right normal-metal electrode by the tunneling. Accordingly, we 
must take account of /t, /l and /n in the coupling terms. The corresponding correction to 
each coupling term has been neglected in the literatures, 6 ' 9 so we obtain it by extending the 
argument in ref. 12. The result is 

PT(X ' £) = 4 e ^s(0^s^ (X ^J mj) (/FD(e + eVm]) " /FD(e " eVhl]) " 2/T(Xinj ' £)) 

+ 6(X ~ Xdct) (fFB(e + eV dct ) - / FD (e - eV dct ) 
-Kdct 



2/ T (x dc t, e) + /n(0, e) + / N (0, -e)) }, (12) 



Ni(e) (5(x-x in j) 



Pl(X ' £) = 4e 2 iVs(0)^ s | R h - — (/fd(6 + eMnj) + /fd( £ - eMnj) - 2/ FD (e) 

- 2/ L (x in j,e)) 

+ 6{X ~ XdCt \ fFD(e + eV dct ) + / FD (e - eV dct ) - 2/ FD (e) 

-Kdet 



2/L(x det , e) + /n(0, e) - / N (0, -e)) }, (13) 



where A^s(0) and As are the density of states and the cross-sectional area of the supercon- 
ductor, respectively, and i?i n j represents the tunnel resistance of the left junction. We have 
assumed that quasiparticles in the left normal-metal electrode are in equilibrium with the 
chemical potential fis — eV- m y The distribution function /n in the right normal-metal elec- 
trode obeys 

d 2 

D N -^f N (y,e) + I N (y,e,{fa}) + P N (y,e) = (14) 
with the diffusion constant Z?n and the coupling term 

- 2^^i£ i,FDie) - /FD(e + eVda) 

- h(0,e) + h{x dct ,e) + /L(a;dct,e)), (15) 
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where -/Vn(O) and are the density of states and the cross-sectional area of the normal-metal 
electrode, respectively. The collision integral is given by 



In{x, e,{f}) =-2 J deV N (e,e')- 



2T 



0/( 



x, el 



cos h 2 (^^ 



)/( 



X, € 



sinh(^) CO sh(^^) coshf^ii 



(16) 



2T 



where a nonlinear term is neglected and on is obtained from the expression of as by the simple 
replacement of as — * a N- Since the behavior of /n near the chemical potential \i — eVdct is 
irrelevant for our calculation of the excess current, we can safely neglect the scattering-in term 
in eq. (16). The collision integral is then reduced to In(x, e, {/}) ~ —f(x, e)/r^(e, V^ ct ), where 
tn(6, Vdet) is the phonon-mediated energy relaxation time defined as 



f . . cos! 

m{e, V det ) = 2 / de a N (e, e ) — , 



cosh(^±|^; 



) cosh( 



-eV A . 



(17) 



2T I 2T 1 

We numerically obtain /t(x, e) at x = x^et an d /n(2/ 5 e) at y = as functions of V^et an d 
Vi n j. The excess current I ex can be obtained by substituting the resulting distributions into 
eq. (6). We assume that fT,h(x, e) = at the right end (x = Lg) of the superconductor while 
dfT,h(x,e)dx = at the left end (x = 0). For the normal-metal electrode, we also assume 
that /]\f(y, e) = at the lower end (y = Ln) and df^(y,e)dy = at the upper end (y = 0). 
We set L s = 120 //m, L N = 22 jum, A s = 0.13 x 0.025 ^/m 2 and v4 N = 0.13 x 0.035 /mi 2 . 
The following parameters are employed: T = 0.1 K, A = 2.55 K, D$ = 5.2 x 10 9 /um 2 /s, 
L> N = 4.1 x 10 9 ^m 2 /s, iV s (0) = 1.45 x 10 6 K^/im" 3 , N N (0) = 0.79 x 10 6 K^/mi" 3 , a s = 
4.0 x 10" 5 K~ 2 , a N = 6.4 x 10~ 5 K~ 2 , i? inj =6kfi and R dct = 7 kft. The parameters for the 
conversion time are chosen as r _1 = 1000 K and ((a p ) 2 ) = 0.04. 
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Fig. 2. The d-dependence of the excess current at Vd c t = for e|Vi n j|/A = 1.2, 1.4, 1.6, 1.8 from 
bottom to top. 



In the following calculation, the left normal-metal electrode for quasiparticle injection is 
fixed at Xi n j = 20 fim and our attention is restricted to the electron injection case of Vi n j < 0. 
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Fig. 3. The Vdet-dependence of the excess current when d = 2 for e|Vi n j|/A = 1.2, 1.4, 1.6, 1.8 
from bottom to top. 

We first consider the excess current when the location of the right normal-metal electrode for 
detection is varied. Figure 2 shows I ex as a function of the separation d = x^ct — ^inj between 
the two junctions for e|Vi n j|/A = 1.2, 1.4, 1.6, 1.8 when the bias voltage at the right junction 
is fixed at e Vdet /A = 0. We observe that the excess current decays nearly exponentially 
as I cx ~ exp(— d/X). Since we shall see that the excess current at Vdet = is determined 
by charge imbalance, A should be identified as the charge-imbalance relaxation length. The 
relaxation length weakly depends on V\ n j since the conversion time depends on quasiparticle 
energy e as seen in eq. (9). We obtain A « 4.2 /im for e|Vi n j|/A = 1.2 and A ps 6.8 fim for 
e|Vi n j|/A = 1.8. This weak Vlnj-dependence reflects the fact that the conversion time becomes 
longer with increasing quasiparticle energy e from the gap edge. These values roughly agree 
with the experimental one of 3.8 /im. 8 We next consider how the excess current depends on 
the bias voltage V^t at the right junction. The Vdet-dependence of I ex when d = 2 /j,m is 
displayed in Fig. 3, where Vdet is restricted to the low-bias regime in which comparison with 
the experiment is possible. The injection voltages are again e|Vi n j|/A = 1.2, 1.4, 1.6, 1.8. We 
observe that I cx increases with increasing Vdet- This behavior is also in qualitative agreement 
with the experimental result shown in Fig. 4(b) of ref. 8. 

We now explain why the Vdet-dependence arises in the excess current. If only the en- 
ergy imbalance is completely neglected by setting e) = 0, the resulting excess current 
becomes nearly independent of Vdet- This means that the Vdet-dependence of I cx is caused 
by the energy imbalance. However, since the energy imbalance has no direct contribution 
to the excess current as can be seen from eq. (6), the Vdet-dependence should be attributed 
to a nonequilibrium quasiparticle distribution created in the normal-metal electrode by the 
energy imbalance. The charge imbalance plays only a miner role in the creation of nonequi- 
librium quasiparticles in the normal-metal electrode because it decays much faster than the 
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Fig. 4. Schematic picture of the phonon-mediated energy relaxation process in the normal-metal 
electrode when (a) Vd ct = and (b) Vd ct > 0. The final state density for tunneled electrons 
(holes) increases (decreases) with increasing the bias voltage Vdct, so that the corresponding re- 
laxation time decreases (increases). Thus, the density of electrons (holes) decreases (increases) 
with increasing Vdct- This variation straightforwardly results in an increase of the excess current. 

energy imbalance due to the presence of the conversion process. Equation (6) indicates that 
quasiparticles in the normal-metal electrode can contribute to the excess current only when 
/n(0, e) + /n(0, — e) / for e > A. Such a distribution can be created by the tunneling of 
quasiparticles in the energy mode in collaboration with phonon-mediated energy relaxation. It 
should be emphasized that if the energy relaxation due to inelastic phonon scattering is absent, 
the resulting quasiparticle distribution inevitably satisfies /n(0, e) + /n(0, — e) = reflecting 
the asymmetric nature of the energy mode, i.e., e) = — /l(x, — e). The Vdet-dependence 
of the energy relaxation time 7>r(e, Vdct) plays an essential role in our argument. Note that the 
Vdet-dependence is determined by the final state density for phonon-emission processes (see 
Fig. 4). At Vdct = 0, the energy relaxation time satisfies tn(c, Vdet) = tn(— e, Vd e t) since the 
final state density for the electron with energy e is same as that for the hole with energy e. 
We thus find /n(0, e) + /n(0, — e) = at Vdct = 0, and the quasiparticles have no contribution 
to lex- This indicates that I C x at the zero bias voltage is purely determined by the charge 
imbalance. However, if Vdct > 0, the final state density for the electron becomes larger than 
that for the hole, so that tn(e, V^et) < tn(— e, Vdet)- This results in /n(0, e) + /n(0, — e) < 0, 
and we observe from eq. (6) that such a quasiparticle distribution makes a positive contribu- 
tion to lex- Similarly, if Vdct < 0, we find TN(e, Vdet) > Tn( — e>Vd e t) and thus quasiparticles 
negatively contribute to I ex . The above argument indicates that I ex increases with increasing 
the bias voltage V^et reflecting the quasiparticle distribution in the normal-metal electrode. 
We arrive at the conclusion that the Vdet-dependence of I C x is caused by the nonequilibrium 
quasiparticle distribution /n which is created by the tunneling of quasiparticles in the energy 
mode in collaboration with the phonon-mediated energy relaxation process. 

In summary, we have considered quasiparticle charge and energy imbalances in a thin 
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superconductor coupled with two normal-metal electrodes via tunnel junctions. Our attention 
is focused on the bias-voltage dependence of the excess tunneling current arising at one junc- 
tion when charge and energy imbalances are created by quasiparticle injection at the other 
junction. We have numerically obtained the excess current I ex as a function of the bias volt- 
age Vdet> taking account of not only charge and energy imbalances in the superconductor but 
also nonequilibrium quasiparticles in the detection normal-metal electrode. We have shown 
that a nonequilibrium quasiparticle distribution created in the normal-metal electrode by the 
energy imbalance can contribute to the excess current although the energy imbalance itself 
has no direct contribution. We have also shown that I cx at the zero bias voltage is purely de- 
termined by the charge imbalance, while the energy imbalance indirectly causes a nontrivial 
Vdet-dependence of 7<jet- The obtained voltage-current characteristics qualitatively agree with 
the experiment by Yagi. 

The author thanks R. Yagi for valuable comments. 
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